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Abstract: Analytical and experimental studies were performed to investigate the effect of gear rim thickness on 
crack propagation life. The FRANC (FRacture ANalysis Code) computer program was used to simulate crack 
propagation. The FRANC program used principles of linear elastic fracture mechanics, finite element modeling, 
and a unique re-meshing scheme to determine crack tip stress distributions, estimate stress intensity factors, and 
model crack propagation. Various fatigue crack growth models were used to estimate crack propagation life based 
on the calculated stress intensity factors. Experimental tests were performed in a gear fatigue rig to validate 
predicted crack propagation results. Test gears were installed with special crack propagation gages in the tooth 
fillet region to measure bending fatigue crack growth. Good correlation between predicted and measured crack 
growth was achieved when the fatigue crack closure concept was introduced into the analysis. As the gear rim 
thickness decreased, the compressive cyclic stress in the gear tooth fillet region increased. This retarded crack 
growth and increased the number of crack propagation cycles to failure. 
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Introduction: A common design goal for gears in helicopter or turboprop power transmissions is reduced weight. 
To help meet this goal, some gear designs use thin rims. Rims that are too thin, however, may lead to bending 
fatigue problems and cracks. The most common methods of gear design and analysis are based on standards 
published by the American Gear Manufacturers Association. Included in the standards are rating formulas for 
gear tooth bending to prevent crack initiation [1]. These standards can include the effect of rim thickness on tooth 
bending fatigue [2]. The standards, however, do not give any indication of the crack propagation path or the 
remaining life once a crack has started. Fracture mechanics has developed into a useful discipline for predicting 
strength and life of cracked structures. 

Ahmad and Loo [3] applied fracture mechanics to gear teeth to illustrate the procedure and estimate crack 
propagation direction. Honda and Conway [4] also applied fracture mechanics to simulate tooth crack propagation, 
compute threshold loads, and calculate tooth life. Flasker and Jezemik [5] applied fracture mechanics to gear 
teeth to estimate stress intensity factors and gear life. Researchers at Tohoku University in Japan performed a 
series of analyses and experiments to determine the effect of residual stress on crack initiation and propagation [6,7]. 
Also, Daniewicz, et al. [8] developed a comprehensive, self-contained analysis package to refine the spur gear 
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bending fatigue theory using fracture mechanics. Lastly, Flasker and Pehan [9] described their method for 
calculating crack propagation in gear teeth using fracture mechanics. Much of the work of the above references 
considered only an initial crack and propagation paths were not considered. Many of the references that did 
consider crack propagation assumed the propagation occurred in a straight path. In addition, experimental validation 
of the cited analyses was sparse. Finally, no work using fracture mechanics was performed for thin-rim gears. 

The objective of this study was to determine the effect of gear rim thickness on crack propagation life. From an 
extensive study [ 10 ], linear elastic fracture mechanics was used to analyze gear tooth bending fatigue in standard 
and thin-rim gears. Finite element computer programs were used to determine stress distributions, estimate stress 
intensity factors, and model crack propagation. Various fatigue crack growth models were used to estimate crack 
propagation life. Experimental tests were performed to validate predicted crack propagation results. 


Fatigue Crack Growth: Many machine elements, such as gear teeth, are cyclicly loaded in application. The 
overall fatigue life of such components may be represented by three distinct phases: 1 ) crack initiation, 2 ) crack 
propagation, and 3) final failure. Once crack initiation has occurred, fracture mechanics may be used to estimate 
crack propagation fatigue growth rate and time to final failure. 


The most universally used method to calculate crack propagation fatigue crack growth was postulated by Paris 
and Erdogan [11], Considered were purely mode I loaded specimens subjected to cyclic load. Further considered 
was unstable crack growth such that the stress intensity factor grew with increasing crack size. Paris postulated 
that the rate of crack growth with respect to number of stress cycles was a logarithmic relationship with the stress 
intensity factor range as , 

-■7 = C(AK) n (1) 

dN 

where da is the change in crack length for dN number of stress cycles, AK is the range of the mode I stress 
intensity factor at a given time, and C and n are material constants. The material constants, C and n , must be 
determined by some experimental means. 


Further research of fatigue crack growth has shown that there exists three important factors not considered in the 
Paris model. First was the effect of load ratio, R , on crack growth ( R = minimum cyclic load / maximum cyclic 
load). Second was the instability of crack growth observed when the stress intensity factor range approached the 
material’s fracture toughness index, Kjc . Third was the presence of a stress intensity threshold factor, AK t h. The 
stress intensity threshold factor is the highest stress intensity factor in which no crack growth would occur. The 
Collipriest crack growth model [12] accounts for these effects where 
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In addressing applications to gears, Inoue, et al. [7] describes fatigue crack growth of gear bending fatigue tests. 
Here, crack growth equations were derived as a function of crack depth through a gear tooth. The expression 
derived for crack growth rate da/dn, as a function of stress intensity range, AK, was 


da 


~dN 

A AK^K,^ 


(l-a^A^-AA^) 


for A K th < A K< A K c 
for AK C < AK < A K IC 

where the parameters Kjc, cc , AKc, AK,h, 77 , and A were all estimated as a function of tooth hardness [7], 
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Crack Propagation Simulation: The analysis of the current study used the FRANC (FRacture ANalysis Code) 
computer program described by Wawrzynek [13]. FRANC is a general purpose finite element code for the static 
analysis of cracked structures. FRANC is designed for two-dimensional problems and is capable of analyzing 
plane strain, plane stress, or axi-symmetric problems. 

Among the variety of capabilities, a unique feature of FRANC is the ability to model a crack in a structure. 
FRANC uses a method called “delete and fill” to accomplish this. To illustrate, the user would first define an 
initial crack by identifying the node of the crack mouth and coordinates of the crack tip. FRANC will then delete 
the elements in the vicinity of the crack tip. FRANC will next insert a rosette of quarter-point, six-node triangular 
elements around the crack tip to model the inverse square-root stress singularity [14,15]. Finally, FRANC will fill 
the remaining area between the rosette and original mesh with conventional six-node triangular elements. The 
user can then run the finite element equation solver to determine nodal displacements, forces, stresses, and strains. 

A further unique feature of FRANC is the automatic crack propagation capability. After an initial crack is inserted 
in a mesh, FRANC models a propagated crack as a number of straight line segments. For each segment, FRANC 
models the crack tip using a rosette of quarter-point elements. FRANC then solves the finite element equations, 
calculates the stress intensity factors, and calculates the crack propagation angle. After the crack propagation 
angle is determined, FRANC then places the new crack tip at the calculated angle and at a user-defined crack 
increment length. The model is then re-meshed using the “delete and fill” method described above. The procedure 
is repeated a specific number of times as specified by the user. In the current study, the stress intensity factors 
were determined from the calculated nodal displacements using the displacement correlation method [16]. The 
method of Erdogan and Sih [17] was used in the current study to determine the crack propagation angle. 

Once the stress intensity factors are determined for each segment, the predicted number of crack propagation 
cycles can be estimated using the fatigue crack growth models. Regardless of the model used, the crack growth 
rates, da/dN, were of the form 

M) (4) 

dN 

where g(AK) is given by Eq. (1) for the Paris relationship, Eq. (2) for the Collipriest relationship, or Eq. (3) for 
Inoue’s method. The predicted number of crack propagation cycles for the i,h crack segment, Nj, was estimated by 


a ; - a 
N ; = — 


g(AK;) 


- + *M 


(5) 


where a x was the crack length of the i t h segment, a^\ was the crack length of the (i-l)th segment, iV/.i was the 
number of cycles of the (i-l),/, segment, and g(AK\) was the average crack growth rate of the i th and (M)m 
segments. Note that a\ was the initial crack length, ^1=0, and i varied from 2 to the total number of segments. 

Gear Finite Element Modeling: Basic gear tooth geometry data was input to a tooth coordinate generation 
computer program. The tooth coordinate generator program used the method of Hefeng, et al. [18] to determine 
the tooth coordinates. The output was tooth coordinate and rim coordinate data which defined a single-tooth 
sector of a gear. This output was used by a commercial available pre- and post-processing finite element analysis 
software package [19]. This package created the finite element mesh of the complete gear. FRANC then used this 
mesh and performed crack propagation simulations. 

Figure 1 shows a sample finite element mesh of an uncracked gear. The tooth geometry used modeled that of the 
test gears of the NASA Lewis Spur Gear Fatigue Rig (described in the following section). The analysis used 
8-node, plane stress, quadrilateral finite elements. The mesh was refined in the region of the loaded tooth for 
improved accuracy. The model of Fig. 1 had 2353 elements and 7295 nodes. Material properties used were that of 
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AISI 9310 steel. The tooth load was placed at the highest point of single tooth contact. For boundary conditions, 
four hub nodes were fixed. In addition, gears with various rim thicknesses were modeled. The parameter describing 
the rim thickness was the backup ratio, ms, where 



where b was the rim thickness, and h was the tooth whole depth. Gears with various backup ratios were modeled 
by incorporating slots in the model. All cases used the same finite element mesh for the loaded tooth. 

Test Facility: Crack propagation experiments were performed in the NASA Lewis Spur Gear Fatigue Rig (Fig. 2). 
The test stand operated on a torque-regenerative principle in which torque was circulated in a loop of test gears 
and slave gears. Oil pressure was supplied to load vanes in one slave gear which displaced the gear with respect 
to its shaft. This produced a torque on the test gears, slave gears, and connecting shafts proportional to the amount 
of applied oil pressure. A 19 kW (25-hp), variable-speed motor provided speed to the drive shaft using a belt and 
pulley. The lubricant used for the gears, bearings, and loading system was a synthetic paraffinic oil. The test gear 
lubricant was filtered through a 5-micron fiberglass filter. 

Test Gears: The test gears were 28-tooth, 8-pitch, 20° pressure angle external spur gears with a face width of 
6.35 mm (0.25 in.). The teeth had involute profiles with linear tip relief starting at the highest point of single tooth 
contact and ending at the tooth tip at an amount of 0.013 mm (0.0005 in.). All test gears used in the experiments 
were fabricated and machined from a single batch of material. The test gear material was consumable-electrode 
vacuum-melted AISI 9310 steel. The gears were case-carburized and ground. The teeth were hardened to a case 
hardness of Rc 61 and a core hardness of Rc 38. The effective case depth (depth at a hardness of R c 50) was 
0.81 mm (0.032 in.). Two different test gear designs were considered. The first was a thick-rimmed gear with a 
backup ratio of m^= 3.3 (Fig. 3(a)). The second was a thin-rimmed gear which incorporated slots (Fig. 3(b)). The 
backup ratio of the thin-rimmed gear was ma=0.3. 

It was believed that tooth bending fatigue cracks would be difficult to initiate based on the load capacity of the test 
rig. Due to this, notches were fabricated in the fillet region (loaded side) on one tooth of each of the test gears to 
promote crack initiation. The notches were fabricated using electrodischarge machining (EDM) with a 0.10-mm 
(0.004-in.) diameter wire electrode. The nominal notch dimensions were 0.20 mm (0.008 in.) in length and 0. 1 3 mm 
(0.005 in.) in width along the full face width of the tooth. The notches were located at the same location for both 
test gears. This location was at a radius of 40.49 mm (1.594 in.) on the fillet, which was the position of the 
greatest tensile stress for the solid gear (m^=3.3). The notches produced a stress concentration factor of 
approximately three as determined using a finite element analysis. 

Instrumentation: The standard test rig instrumentation monitored test gear speed, oil load pressure, test gear 
and slave gear oil pressure, and oil temperatures. Also, overall test stand vibration was monitored using an 
accelerometer mounted on the top housing. In addition to the standard facility vibration sensor, an advanced 
vibration processing diagnostic system was installed in the test stand to help assist in crack detection. Crack 
propagation gages were used in the experiments to determine fatigue crack growth. Special gages were fabricated 
for installation in the tooth fillet region of the test gears. The gages had ten circular strands with an inner radius of 
1 .52 mm (0.060 in.) and an outer radius of 3.05 mm (0. 1 20 in.) (Fig. 4). The strands were designed to break as the 
crack propagated though them, which in turn, increased the electrical resistance of the gage (Fig. 4(a)). Figure 4(b) 
shows the installation of a gage in the fillet region of a notched tooth. A gage was installed on each side of the 
tooth flank for each gear instrumented with crack gages. The electrical resistance of the crack gages were monitored 
along with the load cycle count to estimate cycles as a function of crack length. The information from the rotating 
crack gages was transferred through brush-type slip rings. Also, an infrared tach sensor was used to measure 
number of load cycles. 

Measured Gear Fatigue Crack Growth: The thin-rimmed gear was used in test 1 . The test was run at 89 N-m 
(786 in.-lb) torque and 10,000 rpm speed for 6.5 hr, at which time rim fracture occurred. Figure 5 plots the 
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number of load cycles as a function of the measured crack length. The crack gage results indicated the crack 
growth was non-uniform throughout the tooth face width. A crack started on the rear flank of the tooth at the tip 
of the notch and reached an initial size of 0.46 mm (0.01 8 in.) at 1 ,060,000 cycles. The crack continued to propagate 
through the rear flank but did not reach the front flank until approximately 2,680,000 cycles. At 2,910,000 cycles, 
the crack reached a size of 0.64 mm (0.025 in.) on the front flank, but completed propagated through the rear gage 
by this time. Even though the crack initiation time was not uniform throughout the tooth face width, the crack 
propagation rate was uniform. This was indicated by the similarity in slopes of the curves in Fig. 5 for gages 1 
and 2. 

The thick-rimmed gear was used in test 2. This gear was run at 136 N-m (1200 in.-lb) torque and 10,000 rpm 
speed for 15 minutes, at which time tooth fracture occurred. Figure 6 gives the processed crack propagation 
results for test 2. Note that the crack initiation and crack propagation was fairly uniform throughout the tooth face 
width for this test. 

Comparison of Predicted and Measured Crack Growth: The FRANC computer program was used to simulate 
crack propagation and calculate mode I stress intensity factors as a function of crack length. The predicted stress 
intensity factors were then used with three different fatigue crack growth models (Paris, Collipriest, and Inoue) to 
estimate crack propagation life. 

A comparison of predicted crack propagation cycles using the Paris, Collipriest, and Inoue methods is shown in 
Fig. 7. For this, the thin-rimmed model (ms= 0.3) was used to simulate the test gear of Fig. 3(b). An initial crack 
of 0.64 mm (0.025 in.) was placed in the tooth fillet at the location of the maximum tensile stress. Crack propagation 
was then simulated and the mode I stress intensity factor as a function of crack length is given in Fig. 7(a). From 
this, six different fatigue growth cases were considered. The first four cases used the Paris equation and material 
constants of AISI 9310 specimens from experiments of Au and Ke [20]. The fifth case used the Collipriest equation 
and AISI 9310 material constants from Forman and Hu [21]. The load ratio used was R=- 2.6 as determined from 
the finite element analysis. The sixth case used Inoue’s method and the material constants of the SCM415 material 
(SCM415 is a high-strength Japanese steel, similar in properties to AISI 9310). The predicted number of cycles 
per crack length varied significantly among the cases studied (Fig. 7(b)). Note that the cycles were defined as the 
number of crack propagation cycles after an initial crack of 0.64 mm (0.025 in.). 

Predicted crack growth for the mg= 0.3 and 3.3 gears were compared to the measured crack growth from the 
experiments. Again, the six different prediction schemes as mentioned above were used. The predicted number of 
crack propagation cycles using the sixth schemes were, on the most part, extremely low compared to the measured 
number of cycles from the experiments. To account for this, the concept of fatigue crack closure was investigated. 
Elber [22] performed crack experiments on aluminum alloys and deduced that residual compressive stresses 
existed near the crack tip region due to plastic deformation. These residual stresses reduced the effective stress 
intensity factor range (and thus, increased crack propagation life) and provided a better fit to experimental data 
than other empirical expressions. Elber proposed an effective stress intensity range ratio, U, such that 

AK gJf =U(AK ) (7) 

where AK e jf was the effective stress intensity factor range. Elber then used the effective stress intensity factor 
range in the Paris fatigue crack growth model. In addition, Elber defined U through experimental studies as a 
linear function of the load ratio, R. 

The concept of fatigue crack closure was applied to the current gear crack experiments and predictions. A study 
was then conducted to estimate the effective stress intensity factor range ratio for the gear crack experiments. The 
predicted number of crack propagation cycles using the previously mentions six schemes were plotted versus 
crack length at a variety of arbitrarily chosen U ratios. For the Paris equation and material constants n=2.954 and 
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C=8.433 x 10-9 mm/cyc/(MPaVm)n, good correlation between predicted crack cycles and the experiments occurred 
when: 1) I/=0.4 for R=- 2.6, and 2) C/=0.8 for R=-0. 1 . Assuming a linear relation between U and R produced 

U = 0.82 + 0.]6(R) (8) 

Figure 8 shows a sample comparison of predicted and measured crack growth when the fatigue crack closure 
concept was used. The cycles were defined as the number of crack propagation cycles after an initial crack of 

0.64 mm (0.025 in). It should be noted that good correlation was also achieved when the Collipriest equation was 
used with certain U values. This produced a relationship similar to Eq. 8 but with different coefficients [10]. 

Figure 9 displays the effect of rim thickness on predicted mode I stress intensity factors and predicted crack 
propagation cycles. The stress intensity factors were determined from FRANC using the appropriate finite element 
models. The Paris equation was used along with the effective stress intensity range ratios of Eq. 8. The initial 
cracks of the various models were placed at the location of the maximum tensile stress in the tooth fillet. The 
stress intensity factors were lowest for the ms— 0.5 case. This gave the highest predicted number of cycles for the 
cases studies. The cycles all were defined as the number of crack propagation cycles after an initial crack of 

0. 28. mm (0.01 1 in). The stress intensity factors were highest for the mg= 0.3 case. However, the predicted life for 
this was somewhere between the case of /ng= 0.5 and 1.0 due to the fatigue crack closure effect. The cases of 
rriB= 3.3 and 1 .0 gave nearly the same predicted life. 

Conclusions: Analytical and experimental studies were performed to investigate the effect of gear rim thickness 
on crack propagation life. The following conclusions were made: 1) Good correlation between predicted and 
measured gear crack growth was achieved when the predictions used the Paris crack growth equation and the 
concept of fatigue crack closure. 2) For thin rims, a decrease in rim thickness caused an increase in both the stress 
intensity factor and the compressive cyclic stress in the gear tooth fillet region. The increase in stress intensity 
factor promoted crack growth while the increase in cyclic compressive stress tended to retard crack growth and 
increase the number of propagation cycles to failure. 
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Figure 2. — NASA Lewis Spur Gear Fatigue Rig. Figure 3. — Test gears used to determine effect of rim 

thickness on crack propagation, (a) m B = 3.3. (b) 
m B = 0.3. 
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Figure 5.— Crack propagation fatigue growth for test 
1 , m B = 0.3. 
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Figure 6. — Crack propagation fatigue growth for test 
2, m B = 3.3. 


Figure 4.— Specialized crack propagation gages for 
gear tooth crack growth measurements, (a) In- 
crease in gage electrical resistance as the number 
of broken strands increase, (b) Installation of crack 
propagation gage on test gear. 
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Figure 7. — Comparison of predicted crack propaga- 
tion cycles using Paris equation, Collipriest equa- 
tion, and Inoue equation. Model for m B = 0.3. 

(a) Mode I stress intensity factor, (b) Life 
comparison. 



Figure 8. — Comparison of predicted crack propaga- 
tion cycles to experiments. Paris fatigue crack 
growth model, n = 2.954, C = 8.433 x 10 -9 mm/cyc / 
(MPaVm) n , R = -2.6, U = 0.4, used for predictions. 
Test 1 , gage 1 , front flank, for experiments. 




Crack length, mm 

Figure 9.— Effect of backup ratio on stress intensity 
factors and crack propagation cycles, (a) Mode I 
stress intensity factors, (b) Crack propagation 
cycles, Paris fatigue crack growth model, n = 2.954, 
C = 8.433 x 10 -9 mm/cyc/(MPaVm) n , U = 0.82 + 

0.16 R. 
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